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Abstract 

Comparisons of multiplicity distributions at ^/s = 1.8 TeV (E735 Experiment and 
CDF Collaboration) with our predictions by a two-component stochastic model 
have been made. This stochastic model is described by the pure-birth and Poisson 
processes. It is found that there are discrepancies among data at the CERN SppS 
collider and those at the Tevatron collider. The latter data by E735 do not contain 
the leading particle effect described by the Poisson process, providing that the view 
of the two-component stochastic model is correct. The data by CDF contain small 
leading particle effect. The reason is probably attributed to the correction made by 
E735, which is necessary to make the data of the full phase space from the restricted 
rapidity < 3.25. 
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1 Introduction 



In 1993, we made predictions for multiplicity distributions at the Tevatron 
collider by a two-component stochastic model including the pure-birth (PB) 
process and the Poisson process [1]. See also [2]. Those predictions have been 
based on analyses of Cq = {n'^) / ijiY (g = 3 ~ 5) in the energy range y/s = 
11.5 ~ 900 GeV [3-12]. The method of an extrapolation has been utilized for 
the predictions. 

On the other hand, recently Alexopoulos et al., the E735 Experiment at the 
FNAL [13] has reported the multiplicity distributions of the full phase space at 
300, 546, 1000 and 1800 GeV. The data corrected by the simulation program 
are shown as the data of full phase space. Thus it is possible to compare our 
predictions with the multiplicity distributions at the Tevatron collider. Indeed 
we are very interested in these comparisons, because we can confirm whether 
or not the hadronization process is ruled by the two-component stochastic 
model. 

In the next paragraph, the two-component stochastic model is introduced. In 
the 3rd paragraph, we compare the predictions with observed data by E735 
at the Tevatron at ^/s = 1.8 TeV. In the 4th one, we consider why there are 
discrepancies among the predictions and the data by E735 at the Tevatron. 
In the final one, we present our concluding remarks. 



2 Two-component stochastic model 

We explain bricfiy the essentials of the two-component stochastic model used 
in ref. [1]. It is based on the following two-component branching equation for 
the two-component probability P{na,nh;t): 



In eq. (1) the Poisson component (type a particles) and the pure-birth (PB) 
component (type b particles) do not couple to each other: The former is mainly 
related to particles in the fragmentation region whereas the latter is related 
to those in the central region. The two-component probability P{na,nb;t) is 
therefore a product of multiplicity distributions for the Poisson component, 
Pai'iT'a'jt) and the corresponding one for the PB component, Ph{nb;t): 



dP{na, rib] t) 
dt 



II [P{na - 1, Ub; t) - P{na, Ub; t)] 

+X [{ub - l)P{na, Ub -l;t) - ribPiria, Ub, t)] . 



(1) 



P{na,nb;t) = Pa{na;t)Pb{nb;t), 



(2) 
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and the multiplicity distribution in the final state (i.e., for maximum t — T 
(the maximum value of an evolution time)) is 

P(n)= Yl Pina,n,;T), (3) 

na+ni,=n 




t = 



Fig. 1. Illustration of the two-component stoehastic model: The dashed line denotes 
the virtual source (S) which cannot be observed. The particles a are produced by 
the Poisson process, and b by the PB process. 

Here we use the following initial condition [2]: 

P{na,n,-t^ 0) = ^e-"-^e-"^°, (4) 



which leads to 

Pa(na,(n,)) = ^^e-<->, (5) 



and 



Pb{nb, {ub)) 



{rib) P 



rib 



p-nbo J (1) 




for Ub — Q 
for Ub > 1, 



(6) 



where p = e^^ — 1. Here and after L!^\—x) denotes the normalized Laguerre 
polynomials. The total multiplicity distribution is therefore given by 
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Pi'Ti) = Pain-a, {na))Pbinb, {rib)) 

na+ni,=n 



— g-(("a>+n6o) 



hiin- jy-\'i-+p) j P ^~^\ p) 



(•7) 



The corresponding factorial moments for this P{n) are given by 
F(') = (n(n - 1) ■ ■ ■ (n - / + 1)) 



(8) 



The Cg moments which we calculate in what follows are easily derived from 

F«. We use the same parameters: fiT. n„o- AT and n^o [1] which were deter- 
mined by fits to data for energies y/s = 11.5 ~ 900 GeV [3-12] by choosing 
minimum chi-squared of sum of Cg (g = 2 ~ 5). Finally we have the following 
expressions: 



^ = exp(AT) = exp (x In J^) , (9) 
(na)=//T + n„o, (10) 

with Ubo = 4.2, A = 0.42 ± 0.01, = 8.60 ± 0.33 GeV, /i = 1.22 ± 0.10 and 
UaO = 1.50 ±0.12. 



3 Compcirisons of predictions with data at ^/s — 1.8 TeV 



Using eqs. (9) and (10), wc can predict (n) and Cg = (n^) / {n)"^ in Table 1 and 
Fig. 2. Combining the data at y/s = 1.8 TeV in ref. [13] and our predictions, 
we find that there are discrepancies among our predictions and the data. 

Table 1 

Prediction for multiplicity by eqs. (9) and (10) and corrected data of the full phase 
space at the Tevatron 





in) 


C2 = (n2)/(n)2 


prediction [1] 


47.88 ± 1.54 


1.316 ± 0.007 


corrected data at s/s = l.i 


I [13] 45.81 ± 0.77 


1.45 ± 0.05 



To elucidate the reasons, we adopt other ways. Using eq. (7) with the CERN- 
MINUIT program, and the multiplicity distributions at the FNAL, the CERN 
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n 



Fig. 2. Prediction for multiplicity by eqs. (9) and (10) and data at the Tevatron. 
The solid line is obtained by (n) = 47.88 + 1.54 and C2 = 1.316 + 0.007. The dashed 
line is obtained by (n) = 47.88 - 1.54 and C2 = 1.316 - 0.007. 

SppS and the Tevatron, we obtain various results given in Table 2. See Fig. 3. 

Table 2 shows that the data at the Tevatron do not contain the multiplicity 
by the Poisson process in eq. (7), because (n^) = 0. On the other hand, the 
leading particle effect (finite (ria)) can be observed in the data by UA5 and in 
the ISR regions. 



4 Differences between data by UA5 and those at Tevatron 

To consider the reason of {jig) = at the Tevatron, we analyse the data at 
^/s = 546 GeV and ^=1.8 TeV by CDF Collaboration [14]. What we expect 
is small (na)'s, because of no correction by the simulation program in ref. [14]. 
Indeed we find that the magnitude of (ua) is small, due to the restriction of 
the pseudo-rapidity cutoff \r]\ < 3.25. This fact suggests that the data by CDF 
mainly produced in the central region. (See also [15].) 

As the next step, we carefully compare the data at y/s = 546 GeV by UA5 
and E735. We find that the comparison of the data {^/s ~ 546 GeV) of UA5 
and E735 shows us that error bars of low multiplicities SP{n) {n — 2, 4, 6) at 
the Tevatron are larger than those by UA5. Then we can examine the above 
observation by making a pseudo-multiplicity at ^/s = 546 GeV as follows: 
The magnitude of error bars 5P{n) at n = 2, 4 are assumed to be 5P{n) ~ 
P{n) [n = 2, 4) in Table 4, to consider our problem. As the same situation 
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Table 2 

Analyses of data [3-8,?, 10-13] by eqs. (7) (10) with the CERN MINUIT program. 
(n) = (ua) + {rib) with (nfe) = infeo[l + yjl + 2(n)^(C2 - 1 - l/(n))/nfco]. The data 
by UA5 [10-12] are arranged by means of <tnsd = o"nd + cdd, where NSD, ND and 
DD stand for "non-single diffractive" , "non-diffractive" and "double diffractive" , 
respectively. 



bxps. 








(n) 




C2 






2 / AT 


FNAL 300 GeV 


5, 


.301±0. 


,423 


8.433 


1 


.259±0, 


,028 





24.2/11 


FNAL 800 GeV 


5 


.370ib0. 


,053 


10.16 


1 


.274ib0, 


,004 





48.1/13 


ISR 30.4 GeV 


2, 


.711±0. 


,611 


10.58±0.13 


1 


.186±0, 


,006 


5.259 


19.5/14 


ISR 44.5 GeV 


2, 


.744ib0. 


,541 


12.11ib0.11 


1 


.198±0, 


,007 


5.731 


5.12/16 


ISR 52.6 GeV 


4, 


.366±0. 


,677 


12.73±0.10 


1 


.206±0, 


,005 


3.487 


4.73/18 


ISR 62.2 GeV 


6, 


.540±1. 


,417 


13.67±0.129 


1 


.194±0, 


,005 


1.212 


22.5/17 


UA5 200 GeV 


3, 


.419±0. 


,592 


21.46±0.29 


1 


.249±0, 


,010 


7.005 


9.58/28 


UA5 546 GeV 


3, 


.360±0. 


,180 


29.52±0.18 


1 


.277±0, 


,004 


8.893 


53.0/44 


UA5 900 GeV 


2, 


.703±0. 


,160 


35.93±0.33 


1 


.292±0, 


,007 


13.08 


55.6/51 


E735 300 GeV 


5, 


.963±0. 


,018 


25.97 


1 


.297±0, 


,001 





65.8/57 


E735 546 GeV 


4, 


.993±0, 


,020 


31.06 


1 


.368±0, 


.002 





49.2/77 


E735 1000 GeV 


4, 


.699±0. 


,020 


38.96 


1 


.400±0, 


.002 





112.9/74 


E735 1800 GeV 


4, 


.565±0, 


,002 


46.36 


1 


.417±0, 


.000 





330.6/122 



Table 3 

Analyses of data [14] by eqs. (7) ~ (10). Data are read by eye-balls from their Fig. 
10.^ 



Exps. 


"-60 


(n) 


C2 




xVNdof 


CDF 546 GeV 


6.708±0.348 


36.03±0.21 


1.242±0.004 


1.967 


201.2/43 


CDF 1.8 TeV 


4.754±0.150 


45.22±0.40 


1.358±0.007 


2. MA 


157.0/64 



is seen in P{n) at y/s = 1.8 TeV in Table 5 we make two kinds of pseudo- 
multiplicity distributions and analyse them. 



We have confirmed that the multiplicity distribution derived from the generalized 
gamma distribution in QCD [16] can explain data with smaller chi-squared values 
than those of Table 2. However, it is difficult to find physical meanings of a set of 
parameters: 

00 

P,.(n, (n), D) = ^iV /z"+'^^-ie-^<">-^'*^"dz , 
n! J 


where z = n/{n), N = ^iD^"^ /T{k) and A; = 1/2 {x^/Ndof = 88.8/123 at = 1-8 
TeV). (See also refs. [17,18].) 



6 




50 100 150 200 250 

n 

Fig. 3. (a) Analyses of data at = 546 GeV by eq. (7) with the CERN MI- 
NUIT program, x^/^dof = 49.2/77. (b) Analyses of data at = 1-8 TeV. 
xV^dof = 330.6/122. 

The results of Tables 4 and 5 show that (ua) and C2 depend on the magnitudes 
of error bars at n = 2 and 4, as is shown by analysing the pseudo-P(n)'s. The 
low multiplicities with larger error bars are equal to -P(2) = P(4) = 0, as is 
seen in the first pseudo-P(n) in Table 5. We find that (ua) from the second 
pscudo-P(n) in Table 5 is coincided with the prediction of Table 1. Thus, it 
can be said that the leading particle effect is relating to low multiplicities and 
the magnitude of error bars. 
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Table 4 



P(n) and pseudo-P(ra) with 6P{n) ~ 


P{n) {n = 2, 4) at ^ 


= 546 GcV 




P{n) 


P{n) 


pscudo-P(n) with 




= 546 GeV 


= 546 GeV 


5P{n) - P(n) (n = 2, 4) 




UA5 


E735 Experiment 


= 546 GeV E735 


n 


(real P{n) ± 5P{n)) 


(real P{n) ± 5P{n)) 


(pseudo-P(n) ± 5P{n)) 


2 


0.0027 ± 0.0008 


2.66 X 10-^ ± 0.0997 


2.66 X 10-^ ± 2.66 x 10"^ 


4 


0.0077 ± 0.0013 


0.0201 ± 0.0799 


0.0201 ± 0.0201 


6 


0.0122 ± 0.0014 


0.0423 ± 0.0571 


the same as left column 


8 


0.0195 ±0.0017 


0.0364 ± 0.0496 






53.0/44 


49.2/77 


80.3 


{{ria), {rib)) 


(8.89, 20.63) 


(0, 31.06) 


(8.09, 24.17) 


in) 


29.52 ± 0.18 


31.06 


QO 0(5 -I- n QO 
oZ.ZO it U.oz 




1.277 ±0.004 


1.368 ± 0.002 


1.0U4 ± U.UU ( 


Table 5 








The same as 


Table 4 but = 1-8 TeV 






P(n) 


1st pseudo-P(n) 


2nd pseudo-P(n) with 




\/s = 1.8 TeV 


P(2) = P(4) = 


5P{n) ~ P(n) (n = 2, 4) 




E735 Experiment 


^ = 1.8 TcV E735 


^=1.8 TcV E735 


n 


(real P{n) ± 5P{n)) 


(pseudo-P(n) ± 5P{n)) (pseudo-P(n) ± 5P{n)) 


2 


1.75 X 10"^ ±0.1038 





1.75 X 10"^ ± 1.75 X 10"^ 


4 


0.0223 ± 0.0291 





0.0223 ± 0.0223 


6 


0.0224 ± 0.0327 


the same as left column the same as left column 


8 


0.0167 ±0.0150 








330.6/122 


314.8/120 


432.3/122 


{{ria), {rib)) 


(0, 46.36) 


(0, 46.82) 


(7.94, 39.26) 


in) 


46.36 


46.82 


47.21 ± 0.19 


C2 


1.417 ± 0.000 


1.410 ± 0.000 


1.372 ± 0.003 



5 Concluding remcirks 



We have compared our predictions [1] for the multiplicity distribution at 1.8 
TeV, based on the two-component model, and the data at the Tevatron by 
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E735. It is found that there are discrepancies between them. See Fig. 2 and 
Table 1. 

To elucidate physical reason of the discrepancies, we have analysed the data 
by means of the two-component model. From this procedure, we have known 
that the data by UA5 contain the reading particle effect, finite {na). On the 
other hand, the data by the Tevatron do not contain it [^. 

As the next step, we have resolved the reason why the data by the Tevatron 
do not contain (n^). For this purpose, we analyse the data by CDF [14] and 
find small (ria) at ^/s = 546 GeV and 1.8 TeV. See Table 3. Moreover, we 
compare the data by UA5 and at the Tevatron in Tables 4 and 5 and find 
that the error bars of low multiplicities, 6P{n) {n = 2, 4, 6) are large. To 
see how our result depends on large 5P{n) {n = 2, 4, 6), we have made the 
pseudo- multiplicity distribution, assuming 6P{n) ^ P{n) aX n = 2 and 4. 
From the pseudo-multiplicity distribution, we can estimate the finite leading 
particle effect at -^i = 546 GeV and 1.8 TeV at the Tevatron. The result of 
(ua) ~ 8 from the second pseudo-P(n) in Table 5 is almost the same of our 
prediction [1]. (See Table 1.) 

In conclusion, we have found that the corrected data of full phase space at 
Tevatron by E735 do not contain the leading particle effect which is observed 
in the data by UA5. It is difficult to estimate this effect, because of larger error 
bars at low multiplicities, P{n) ± 6P{n) (n = 2, 4) In other words, the 
low multiplicities with 6P{n) > P{n) {n = 2, 4) are equivalent to P{n) = 
(n = 2, 4). They are playing important roles for determining the leading 
particle effect which is expected in ref. [1]. 

Moreover, providing P{n, 61]) with the pseudo-rapidity cutoff in a future, we 
could obtain more useful information for the production mechanism which are 
illustrated in Fig. 4: The pure-birth (PB) process mainly describes P{n) in 
the central region. On the other hand, the Poisson process mainly does P{n) 
in the fragmentation region [^. 



^ It should be noticed that the streamer chamber at the CERN SppS and the track- 
ing detector at the Tevatron are used. The pseudo-rapidity range of the Tevatron 
is \r]\ < 3.25, which corresponds to the central region. As events of low multiplicity 
probably contain particles with larger r], it may be said that the leading particle 
might be missed in measurements at the Tevatron. Of course, the corrected data of 
the full phase space should contain the effect. 

^ For us, it is difficult to discuss the valuation of the correction methods applied 
to the data at the Tevatron by E735. 

^ Two papers relating to da/dr] are presented in different points of view by Takagi 
and Tsukamoto [19], and Ohsawa [20]. 
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dc 
" dr\ 




Fig. 4. da/dr] distribution. The PB and the Poisson processes probably rule the 
central and fragmentation regions, respectively. 
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